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p Abstract. Single and double inclusive cross-sections in nucleus-nucleus collisions are 

^ ! derived in the perturbative QCD with interacting BFKL pomerons in the quasi-classical 

p . 

approximation. 

m ; 1 Introduction 

a^ . 

^ . With the advent of colliders the study of particle production in nucleus-nucleus collisions 

^ J ! has acquired a prominent role both in experimental and theoretical studies. Single and 

O I double inclusive cross-sections and correlations related to them draw naturally most at- 

tention. In the framework of the perturbative QCD they can be studied in the approach 
based on interacting BFKL pomerons developed for nucleus- nucleus collisions in [H Ej [3] . 
Intuitive considerations plus experience with the Local Reggeon Field Theory (LRFT) al- 
low to guess the structure of the single inclusive cross-section as emission either from the 
I central pomeron or from the two adjacent triple pomeron vertexes in the convolution of 

two sets of fan diagrams propagating from the center to the projectile or target. However 
no rigorous demonstration of this structure has been given in the literature. Still worse is 
the situation with the double inclusive cross-section for which even in the LRFT one has 
a very complicated formula (see [1]). 

In this paper we aim at filling this gap. We derive formulas for both single and double 
inclusive cross-sections in nucleus-nucleus collisions in the perturbative QCD approach 
with interacting BFKL pomerons. Part of the problem which involves particle emission 
from the pomerons is ideologically rather similar to the LRFT (although considerably 
more complicated technically). So to treat this problem we shall use the cut pomeron 
formalism developed within the LRFT [5], which we appropriately generalize for the 
BFKL pomerons. A new problem is emission from the triple pomeron vertex. For the 
single inclusive we find that it is indeed described by simple formulas used in our previous 
studies ( [6]). However the double inclusive cross-sections with a single emission from 
the vertex lead to a complicated expression including different parts of the cut emission 
vertex found in 
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Numerical calculations of the single inclusive cross-section in nucleus-nucleus collisions, 
with vertex emission taken into account, were earlier reported in [6] (as mentioned, 
without rigorous justification of the formulas). As to the double inclusive cross-sections 
and correlations, due to their complexity, we do not attempt here to use them for numerical 
studies. These studies present a separate difficult problem and are postponed for future 
publications. 

The paper is organized as follows. In the next section we generalize the cut pomeron 
formalism for BFKL pomerons. The derivation closely follows that in the LRFT ( [HE]) 
with inevitable complications due to the non-local structure of the BFKL pomeron. Next 
in Section 3 we solve the equation of motion and find working formulas for the single and 
double inclusive cross-section corresponding to emission from the pomeron. In Section 
4 we discuss emission from the vertex and derive formulas for the single and double 
inclusive cross-sections which involve such emission. Finally we discuss our results in the 
last section. 



2 Cut pomeron formalism and inclusive cross-sections 
for nucleus-nucleus scattering in the perturbative 
QCD 

2.1 Fields and Lagrangian 

Pomerons are described by two fields and which depend on rapidity y, relative gluon 
momentum q and point b in the transverse space. We shall restrict ourselves with the 
purely forward case relevant for collisions off large nuclei. Then b is conserved during the 
collision. Correspondingly we do not indicate it explicitly unless it may lead to confusion. 
After transition to real fields the generating functional for the Green functions is written 
as 

' DipDip^e-^, (1) 



where action A is an integral over y, q and b: 

-A = J dyd%d^qL{^,ij^) (2) 



and the Lagrangian density L is 
L(y,g,6)=^tg^ + A[(i 
Here Q is essentially the BFKL operator 



L{y,g,t,) = '*tQV. + \\(IH>I) ■ + (/«.) • '(//'I + gS{y)il, + fS(Y - t,)«)', (3) 



with H the BFKL Hamiltonian, and 

K = Vyvl (5) 

Operator K is conformal invariant and commutes with H. The triple-pomeron coupling 
A is 

A = (6) 

TT 
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Finally (?(?/, b) and f{y, b, B) represent the coupling of the pomeron to the colliding nuclei. 
giy, b) describes this coupling to the target A (at zero rapidity) 

g{y,b)=ATA{b)p{q), (7) 

where T^(&) is the profile function of nucleus A with its center at the origin in the trans- 
verse plane and p(g) is the colour density of the nucleon. /(y, 6, B) describes the coupling 
to the projectile at rapidity Y and its center at point B (the impact parameter of the 
collision) 

f{yAB) = BTB{B-b)p{q). (8) 

To study different cuts we introduce the cut pomeron formalism following the scheme 
developed in [5] for the local pomeron. We introduce 6 fields 

^„ 4, (9) 

which describe pomerons to the left (right) of the cut (subindex +(-)) and cut pomerons 
(sub index c). To fulfill the AGK rules the Lagrangian of these fields is to be taken as 

L, = Ll + Li + Ll (10) 

Here the unperturbed Lagrangian is 

= E ^^4Qi'^. (11) 
i=±,c 

with 

e± = -ec = 1. 
The interaction involves several terms 

Li = X{Kijl) ■ i^l + X{Kij+) ■ i^l" + A(ir^l) ■ ^! + A(ir^_) ■ 

-2\{K^l) ■ + - 2X{K^,) ■ + ^l) 

+ V2\{Kijl) ■ iij+i:. + ijl) + V2X{Kij,) ■ iijlijl + i:f). (12) 
Finally the external Lagrangian is 

Lf = g5iy)i^+ + ^_ - v^^,) + /5(r - y){^l + ^1 - v^^l). (13) 
2.2 Inclusive cross-sections 

To study inclusive cross-sections we additionally introduce interaction of the cut pomeron 
with the emitted particles. As is well-known, the inclusive cross-section corresponding to 
emission of a gluon jet from the pomeron can be written as a double integral 

Hy, 1^) = ^'^^y^2^ = 2 y d%d^qi£q25^{qi + q2 - fi)v{Qi, q2)PA{y, Qi, b)pB{y, q2, b), (14) 
where PA{B){y,q,b) is the pomeron coupled to the target (projectile) and 

V{(lu(l2) = giVi^aVa- (15) 
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This expression can be conveniently rewritten in the coordinate space in the local form 
in terms of the 'non-amputated' pomeron 



P{y, r, b) = 27Tr' J ^e^'^>(y, q, b) 



as 



/ ,21,^2 



d'bd'rPAiy, r, 6)V'e*^'^V'PB(?/, r, b). 



(16) 
(17) 



To generate inclusive cross-sections we add a term to the action 



-AAr 



dyd bd qid q2^{y, qi, q2,b)ipl{y, qi,b)^/jciy, q2,b). 



Differentiation in ^ then gives insertions into the cut pomeron propagator at rapidity y 
and transverse point b with momenta qi and q2- As in [4J we denote the total action as 



AciO = Ac + AAc. 

It follows from unitarity that 

Ac{i = 0) = 

Also one finds that the generating functional of the amplitudes at ^ 7^ is given by 

T,(0 = l-5 + -(e"^^(«)-f 



(19) 
(20) 

(21) 



where S is just the S matrix at ^ = [4J (also see Appendix 1). 
The single and double inclusive cross-sections are obtained as 



hijj, k) = 2 y" d%d^qid^q25'^{qi + q2 - K,)rii(li, ^2] 



S^{y,(lu(l2,b) 



(22) 



and 



hiy, K,\y', k') = 2 J dVb'd'^qid'^q26'^{qi + q2 - K,)d'^q[d'^q'26'^{q[ + q2 - K,')r]ic[i, (l2)viQi, 0.2) 



5ay,qi,Q2,b)day',qW2,b') ?=o- ^^^^ 

Next we use the property proven in [3] that at any values of C,, due to the equations 
of motion, 



■^l{y,(ihb,0My,(i2,b,0- 



(24) 



S^iy,qi,q2,b) 

This together with (120|) immediately gives a simple formula for the single inclusive cross- 
section: 

h{y, = / dVqid^q25'^{qi + q2 - q)viQi, 92)^l(z/, qi,b,^ = 0)i)c{y, q2, b,^ = 0). (25) 
The formula for the calculation of the double functional derivative at = becomes: 



S^{y,qhq2,b)6^{y',q{q2,b' 
+ 



1 r 



?=o 2 
6 



i^l{y, qi, b, OMy^ 92, 0-^liy'^ q'l^ b\ OMy'^ q2^ b\ Ci 



S^{y,qi,q2,b) 



l{y',q[,b',0My',q2,b',0 



(26) 



So the whole problem is reduced to finding ipc, 'ipl and their first derivatives in ^ at ^ = 0. 
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2.3 Unitary transformation to new variables and equations of 
motion 

Considerable simplification can be achieved if, following [4], we make a unitary transfor- 
mation of our fields introducing new fields 0±, 0o and conjugate fields n±, ttq by 

1 ^ 

^± = 0±, ^/'c = --^(00 - 0+ - 0-), 

4 = 7r± + 7ro, = V2no. (27) 
This transforms Lc in to a new Lagrangian 

i=+,-,0 i=+-,0 

+ X{7T+7To ■ K{(j)o + 0+ - 0-) + vr_7ro ■ K{(j)o - 0+ + 0-)) + 7r+7r_ ■ /C(0+ + 0_ - 0o)) . (28) 
The coupling to the nuclei takes the form 

L'e = 95iy)<Po + my - Y)in+ + 7r_) (29) 
and the part of action depending on ^ becomes 

-AAc = y" (i?/rf^gic/^g2^(l/,gi,g2)vro(gi)(0+(g2) + 0~(g2) - 0o('?2))- (30) 

Variation of the action with respect to our fields gives equations of motion. It can be 
seen that they are totally identical for ± fields. So it is sufficient to write them down for 
= 0+ = 0_, 00, vr = 7r_|_ = 7r_ and ttq They are (at point y, q): 

ttQ + 2XK{7i^)) + 2\<pK^ + j (fk^iy, q, k)7io{k) = 0, (31) 

TToQ + XK{7vl - tt" + 27ro7r) + 2A0o ■ - j d^k^{y, q, k)TTo{k) + g6{y) = 0, (32) 
g0 + AK02 + Avr ■ Km - 0o) + Avro ■ if 0o + /5(r - y) = 0, (33) 
Q0O + XK<Pl + 2A(7ro + tt) ■ if 00 + 1 ^2^:^(2/, g, fc) (20(A;) - 0o(fc)) . (34) 

3 Inclusive cross-sections due to emission from the 
pomeron 

3.1 Solution of the equation of motion: fields at ^ = and single 
inclusive cross-sections 

Our first task is to find the solutions of our equations of motion at ^ = 0, which knowledge 
is sufficient for the single inclusive cross-section according to ( l25l) . 

If ^ = then the equations for vr and 0o do not contain driving terms, so that these 
fields are identically zero: 

7r(l/,g) = 0o(l/,g) = 0. (35) 
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The two remaining equations decouple: 

7rog + Ai^(7ro') + (7(5(1/) = (36) 

and 

Qct) + \K(t)^ + f5{Y - y) = 0. (37) 

These are the standard BK equations for the sum of fan diagrams % which go to the point 
(y, q) from the target 

Mv^q) = x{y,(i,9) (38) 

or from the projectile 

Hy,q) = x{Y-y,qJ) = x- (39) 

Put into the expression fl2S]) they give the commonly used factorized expression for 
the single inclusive cross-section as a convolution of two sets of fan diagrams: 

h{y, = 2 y" dVqid'^q25'^{qi + q2 - Q2)x{y, Qi, b)xiy, ^2, b) (40) 

or, in the coordinate space, in terms of Z{y, r, b) related to xiV: b) similarly to f|T6|) : 

72 

Z(y, r, h) = 2nr' J ^e"^''x{y, q, b), (41) 

as 

I{y, «:) = / dVrZ{y, r, b)V^e'^'^V'Z{y, r, b). (42) 

3.2 Derivative of the fields in ^ at ^ = and double inclusive 
cross-sections 

To set up equations for the derivatives of the fields in ^ at ,^ = we have to differ- 
entiate Eqs. fl?Il) - flM|) in ^{yi,qi,q2,bi) and then put ,^ = 0. The field derivatives 
and the equations as a whole will depend on 7 variables: y,q,b and yi,qi,q2,bi, which 
we shall show explicitly only when it is necessary. Thus say Stt/6C, will in fact mean 
Sn{y, q, b)/6^{yi, qi, q2, &i). With these notations we get the equations: 

(- — + // + Ax)2ir^ + 6{y - y,)5\b - b,)6\q - q,)x{q2) = 0, (43) 



2i^(-^ + i/ + Ax)^ + 2A^irx + 2Air(x^)-%-y05'(&-&i)^'(g-giMg2) = 0, 
^ (44) 

{^ + H + Ax)2ir^ + +25{y - y,)5\b - b,)5\q - q,)xiq2) = 0. (46) 

To simplify these multivariable equations we first note that the b dependence of the 
derivatives is trivial: obviously they all are proportional to S'^{b — bi). So we separate this 
factor and consider the derivatives at a fixed point b which need not be shown explicitly. 
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Next we return to our expression for the double inclusive cross- sect ion. Obviously it 
consists of two terms. One is just the product of two single inclusive cross-sections 



The other term comes from the second term in ( 126|) and has the form 
l2\y, f^\y', f^') = J (fh(fh'(fqi(fq25'^{qi + q2 - K)(fq[d'^q26'^{q[ + q'^ - K)r]{qi, q2)ri{q[, ^2) 

5 



(47) 



:i^l{y\q'i.h\i)Uy\q'2.b'i) 



S^{y,(ii,(i2,b) 

J (fqi£q25'^{qi + ^2 - K,)d'^q[d'^q'26'^{q[ ^q^- K')r]{qi, q2)ri{q[, ^2) 



(V^l(i/',gi,&',e = o) 



«=o S^{y,quq2b) 



(4J 



Due to factor 5^(6 — b') contained in the derivatives the double integrations in b and b' 
in fact turn into one over the common point b. Taking this into account and suppressing 
this common argument b, in terms of fields and vr we find 

l2\y, t^\y\ 1^') = j d%d^qid'^q25'^{qi + q2 - K)d^ q^d^ ((2^'^ {q^ + ^2 - 'l2)vic[i, 0.2) 



{x{y\ q'l 



, ,,5(20(y',g^,O-0o(2/',g2,O) 



«=o 5i[y,qi,q2) 



(49) 



5i{.yAiA2) 

where it is assumed that factor 6'^{b — b') has been dropped from the derivatives. From this 
formula we can conclude that we do not need our derivatives in C,{yi, gi, ^2) at all values 
of its arguments but rather integrated over qi and q2 with weight 6{qi + q2 — K^ililii 12) ■ 
Correspondingly we define (suppressing the argument b) 



^o{y,q,yi,ti) 
n(i/,g,i/i,K) 

*o(2/,g,yi,ft:) 



d^qid^q^5^{qi + q2 - K)r]{qi, q2 

d^qid^q^5'^{qi + q2 - /t)?7(gi, ^2^ 
d^qid'^q'^6'^{qi + q2 - n)r]{(li, Q2) 
d^qid^q^5'^{qi + q2 - i^)v{Qi, Q2) 



SMy,Q,0 



S^iyi,qi,q2 

'^7r(?/,g,0 



s^{yi,(ii,(i2) 



H{y,(i,0 



s^{yi,(ii,Q2) 



$(y,9,i/i,/«) = 

In view of expression fl49p we also introduce 

^(y, Q, yi, /«) = 2$(y, g, yi, k) - ^o{y, q, yi, k.)- 

In terms of these quantities we find 

lPiy,^\y',n') 



5=0 



(50) 

(51) 
(52) 
(53) 

(54) 
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(55) 

Integrating our equations fH51) - fHB]) over gi and g2 with weight 5^(gi + ^2 — '^)v{li^ I2) 
we get the following system 

[-^ + H + Xx)2KU + S{y - y^q, k - q)x{y, n - q) = 0, (56) 
2K{--^ + H + Xx)llo + 2X%-Kx + 2XK{xIl)-6{y-yMq,K-q)x{y,K-q)=0, (57) 



dy 

2K[-^ + H + Xx)^ + 8XU ■ Kx + 2XK{x%) - 26{y - yMq, k - q)x{y, K-q)=0, (58) 

{-^ + H + Ax)2ir$o + 25iy - y^q, k - q)xiy, K-q) = 0. (59) 

One observes that the first and last equations determine 11 and $0 in terms of the 
known x and after which the second and third equation allow to find IIo and \Ef. 

Eqs. (!56|) - (!59|) are convenient for numerical calculations. However one can also 
formally express the solution in terms of Green function of the operators d/dy + H + Xx 
and d/dy + H + X ^'^ obtain formulas which can be compared with [4j 

3.3 Formal solution 

We define the Green functions: 



Each Green function is an integral operator in {y, q) space, so that e.g. for G the kernel 
is G{y, q\y' , q'). In our formulas rapidities and momenta enter asymmetrically, so in many 
cases we shall suppress the momenta but leave the rapidities, considering Gyyi as an 
operator in the momentum space with the kernel Gyyi{q\q'). 

In these notation solution of Eqs. fl56l) and fl59l) is immediaite 

Ii{y) = -\K-^Gly^{rix)y., (61) 

%{y) = -K-^Gy,yMx)y.. (62) 
We put these solutions into Eq. (IFTl) multiplied by {2K)~^ to obtain 

(- — + H + Xx)Uo - XK-\Kx{y))K-'Gy,y,{vx)y, 
- \x ■ K-'GUvx)y, - Siy - y,)^K-\vXyJ = 0. (63) 
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(2) 

Figure 1: Diagrams corresponding to the non-trivial part I2 of the double emission from 
pomerons. Crosses mark x ^^^d G, full circles mark x £^nd G. 



G^K-\Kx)K-^G] {r^x)y, + \\[G^xK~'G^ 



yyi 



y,yi 



yi 



GlK-\vx)y.. 



(64) 



Applying operator G^ we get 

no = A 

Similar operations with Eq fl58|) first give the equation 

{^+H+Xx)^-2XK-\Kx)K-'Gly^ivx)y->^xK-'Gy,^^^^^ 

^ . (65) 

which after application of operator G gives 



GK~\Kx) ■ K-'G^ {r,x)y, + A GxK'^G {r,x)y, 
^y,yi L ^y,yi 



^ = 2A 

One observes that is obtained from 2110 by the substitutions 

G^G^, G^G^, x^X- 



Gyy,K r]Xy,. (66) 



(67) 



Putting the obtained expressions for the field derivatives ( l6Tl) . ( l6TI) . ( l64l) and ( l66l) into 
55|) and using the property (1671) we obtain the part 12'^ of the inclusive cross-section as 



li'\y,K\y\K') = 
dVq[d%6'{q[ + q'2- k'Mq',, q'2)Xy',,'^{2X[G^ K-\Kx)K-'G 



y',Q2\y>'i2 



ivx) 



ym 



+ \[g^xK-'G^].,,„ {rix)y,,.+Gl,,,^\y,,,K-\r^x)y,,.] + [G ^ G^ , G^ G^, x x)- 

(68) 

The obtained expression (!68l) is similar to the one in the framework of LRFT [4j, except 
that in our non-local case all quantities are operators also in the momentum space and that 
in various places there appear operators K and acting in this space. Its diagrammatic 
illustration is presented in Fig. [1] In this figure external lines with crosses (circles) show 
sums of fan diagrams propagating towards the target x (projectile x)- Correspondingly 
the internal lines marked with crosses (circles) show the Green functions G (G). Horizontal 
lines indicate the two observed particles. To the diagrams shown in Fig.l one should add 
similar ones with the target and projectile interchanged. 

From the practical point of view expression fl68|) is not very useful due to multiple 
integrations in {y, q) space. Numerical solution of evolution equations (l56l) - (!59l) seems 
more promising. 
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4 Emission from the vertex 

To include emission from the triple-pomeron vertex we have to add new parts to the 
Lagrangian which describe this emission. In the sphtting vertex the pomeron before the 
spht has always to be cut. The two emerging pomerons may be both cut and uncut. In 
accordance with the structure of the interaction in Eq. ( fT2l) the vertex emission part of 
the action is to be taken as 

— Ay = J dyd%(f qi(f q2d'^ q^Ly + /i.e., (69) 

where 

Ly = V2iijl{qi){-fd{qi\q2,q3)M(i2)i^-{(i3) - ^27(91^2, g3)^c(g2)(V^+(g3) +'ip-{q3)) 

+ lc{qi\q2,q3)Mq2)Mq3))- (70) 

The part explicitly shown corresponds to the emission from the splitting vertex. The part 
corresponding to the emission from the merging vertex is indicated as h.c. The vertex 
functions 7^, 7 and 7c describe emissions from the diffractive, single and double cuts 
of the triple pomeron vertex respectively. They are real functions of the three relative 
momenta of the joining pomerons and carry factor + qi — q2 — qs) where k is the 
momentum of the emitted jet. They are different and their form has been found in 
[1] and reproduced in Appendix 2. The fields in fITOl) are supposed to be taken at the 
same rapidity y and transverse point b, which dependence is not shown explicitly. In 
the perturbation expansion the new interaction term has to be taken the number 
of times which corresponds to the number of the emissions from the vertex. So for the 
single inclusive cross-section we have to take it only once and for the double inclusive 
cross-section at most twice. 

The single inclusive cross-section corresponding to the emission from the vertex is 
obtained from (ITOll just by substituting the fields by solutions of the equations of motion 
(pTll - (IMl) at ^ = and integrating over b. For the fields ip and ip"^ these solutions are 



^± = X, = V2x, = X, ^1 = V2X- (71) 
Thus we find a contribution 

Ii^\y,i^) = 2 / d^bY[d^qi(yd{qi\q2,q3) - 47(gi|g2, gs) + 27c(gi|g2, gs)) 
i=i 

{x{qi)x{q2)x{q3) + x{qi)x{q2)x{q3)) ■ (72) 

The total vertex 

7*°* = 7, - 47 + 27c (73) 

has a simple form in the coordinate space. If we introduce the coordinate vertex T{ri \r2, r^) 
acting on non-amputated pomerons according to the relation 



d'^qid^q2d'^q3Xiqihiqi\q2,q3)xiq2)xiq3) = / d^rid'^r2d^r^Z{ri)T{ri\r2,r^)Z{r2)Z{r^), 

(74) 
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where xil) and Z{r) are related by ( HTl) . then one finds [7j 



r*°*(ri|r2,r2) = -^^Vle'^'^'VlS\r2 - n)6\r, - n). (75) 



It corresponds to the expression first obtained in [S] as a contribution additional to the 
emission from the pomeron. 

Passing to the double inclusive cross-section we first find a contribution corresponding 
to the emission from two verteces, which is obtained by taking a product of two interac- 
tions ( |70l) with different external momenta k and k', at different rapidities y and y' and 
transverse points b and b', substituting in them the fields by the solution of the equations 
of motion and integrating over both transverse points b and b'. As a result we obviously 
find a product of two single inclusive cross-sections fl72l) : 



li^'\y, ^\y', k') = li'\y, ^)ll'\y', k'). (76) 

To find the mixed contribution in which one jet is emitted from the vertex and the 
other from the pomeron we have to consider the theory with the interaction term AA^, 
Eq. f|T8|) . and once differentiate in ^: 

I^^\y, K\y'K') = f dVb' n d\d%d%6\q[ + q', - K')i^{q[.q',) 

i = l 

— — -L^(?/, 6, K, qi, ga, gs, Oe"^"^^^] , + (?/ ^ k'). (77) 

Differentiation in ^ will give two terms. One comes from the differentiation of the 
exponential. After integrations over b' and q'l and ^2 with weight 7] this differentiation 
will give the single inclusive cross-section Ii{y',K,') corresponding to emission from the 
pomeron. Factor at ^ = will generate the single inclusive cross-section fl72l) . As a 
result this part gives a contribution 

li'\y, K)h{y\ k') + lt'\y', K')h{y, k). (78) 

If we introduce the total single inclusive cross-section 

lf{y,K)=h{y,K) + I^^\y,K) (79) 

then collecting ( H7I) . (I78il and (1761) we find a factorized contribution to the double inclusive 
cross-section 

lt\y, «:|yV) = P°\y, n)P'^\y'K'). (80) 



The second part of 12^'' will come from the differentiation in ^ of the fields inside L 



7' 



the exponential factor giving unity. Differentiation in ^ together with integrations over 
q[ and q2 with weight r] will substitute fields in accordance with Eqs. (150|) - (154|) . If we 
define 

then in terms of fields ip and we find 



D = f d\d^q2S\qi + q2- /t)r/(gi, 52)777 ^ ^ (81) 



DV',= -^^, D^± = <l>, D^l = V2Uo, Z^V± = n + no. (82) 
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Figure 2: Diagrams for the non-factorized part l'^''^^^'^ of the double emission from a 
pomeron and a vertex. Notations are as in Fig. [1] 

We recall that the derivatives are proportional to 5^(6 — 6'), so that the double integration 
over h and h' turns into a single one. We then obtain the following expression for the 
non-factorized part of the double inclusive cross-section I^^ 



i=l 



+2x(?i)x(?3) 



toil 
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'(?ik2,g3)^(g2) + {ld{qi\q2,q?>) - 27(gi|g2, g3))*o(g2) 

+7*°*(giki,93)^(gi)x(g2)xfe) 
+4x(9i)x(93)[7*°*(gik2,g3)no(g2) + [id{qi\q2^qz) - 27(gi|g2,g3))n(g2)]} 



(83) 



In this expression it is assumed that the derivative fields $0) n, ^ and IIo, apart from 
the argument explicitly shown, depend on their 'own' rapidity y and transverse point h 
and also on rapidity y' and external momentum k' which enter in their definitions (150|) - 
Due to property (1671) the last two terms in (183|) are obtained from the first two ones 
by interchanging the target and projectile. Graphical illustration of l'^^''^^'^ is presented 
in Fig. [2], which shows diagrams corresponding to the six terms in fl83l) . The notations 
are as in Fig. [TJ To the diagrams shown in Fig. [2] one has to add the diagrams which 
are obtained by the interchange of the target and projectile corresponding to the last two 
terms in (l83l) and also diagrams from the the interchange y ^ y', k ^ k' . 



5 Conclusions 

We have derived expressions for the single and double inclusive cross-sections in nucleus- 
nucleus collisions in the framework of the perturbative QCD with interacting BFKL 
pomerons, in the quasi-classical approximation (without loops). The cross-sections in- 
clude terms with emissions both from the pomerons and from the triple pomeron vertex. 
The obtained single inclusive cross-sections are simple. As expected from the AGK rules 
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they reduce to emissions from the central pomeron in the convolution of two sets of fan 
diagrams connecting it with the projectile and target and from the two neighboring verte- 
ces by which this pomeron splits into fans. This expected form of the single inclusive 
cross-sections was used in the calculations in [6]. 

In contrast, the expressions for the double inclusive cross-sections are much more 
complicated. A part of it is just the product of two single inclusive cross-sections corre- 
sponding to independent emissions of two jets from different points in the nuclear overlap 
transverse space. The other part however includes all sorts of rescattering corrections and 
needs summation of diagrams of a structure different from fans. It also includes emissions 
from the triple pomeron vertex having different forms for different cuts passing through 
the vertex. We have set up evolution equations which allow to find the rescattering cor- 
rections to the double inclusive cross-sections. Different forms of emission from the vertex 
have been found in our paper [7]. So in principle the theoretical basis for the calcula- 
tion of the double inclusive cross-section is completed. However practical computations 
along these lines seem to be quite complicated and apparently present a separate (and 
formidable) task. 
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7 Appendix 1. The unitarity equation with cut 
pomerons 

The derivation with BFKL pomerons follows that in [1] for local pomerons. The starting 
point is the equation which tells that the sum of all Green functions in the cut pomeron 
theory is equal to the sum of all Green function in the original theory, provided the 
couplings of the fields to themselves and external sources are chosen in accordance with 
the AGK rules. It is written as 

9_Q(n,m) _Q{ni,mj) 

ml n+!n_!nc! m+!m_!me! ' 

Here G'-"'™'^ are the Green function in the original theory with n incoming and m outgoing 
pomerons and G^""^'"^^^ are the Green function in the cut pomeron theory, where i,j = 
+, — , c. Hi is the number of incoming pomerons of sort i, n = n+ + + Uc and rrij is 
the number of outgoing pomerons of sort j, m = m+ + m_ + rric. Symbol J2' means that 
one should drop terms with n± = Uc = and m± = rric = which are absent in the cut 
theory. The Green functions are assumed to be operators acting in the rapidity-transverse 
momentum space, convoluted with the external sources given by functions g and /. 
One then uses obvious properties 

Q(n±=n,n^=nc=0;m,^=m,m±=mc=0) _ g (86) 
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to consider the sum ( 184|) without restrictions on the summations over Ui and rrij and 



multiphed by 2. One has 

= 2 V ^ ^^"(^/2r-^ ^K,^,) ^/'"(^/2r-'^ ^ ^G^'^^^ + l. (87) 
n.t^onf;;;^, n+!n_!nj m+!m_!mj „,t^o ^' 

where in the first sum of the second equahty n = + + ric and m = m+ + m_ + mc. 
However according to (18^ the first term is exactly equal to the second with the opposite 
sign. So we are left with 

^'."(V2)°;G,„..„,,'r(v^r- 1. (88) 

The sum in (18 BP is just the total S matrix in the cut pomeron theory, so that the equality 
(1881) means that action A^^ = 0) = 0. 

In the theory with the action Ac{C) the amplitude Tc(^) is 

Adding and subtracting terms with n+ = ric = and m_ = vrif. = and using properties 
(IHHD and (1861) we find 



'^^^^ 2 n+!n_!nj ' ^''> m+\mAm,\ ^ n! m! 2' ^ ^ 



Li ^1 I L J 



where again in the first sum n = + n_ + rtf. and m = m+ + m_ + rric- The first sum 
is just one half of e~'^'=*^^^ while the second is e~"^ in the original theory, so that Eq. (190|) 
leads to Eq. 



8 Appendix 2. Cut vertex functions 7^/, 7 and 7c 



The three cut vertex functions 7^, 7 and 7c which enter (jTOj) can be most conveniently writ- 
ten in the coordinate space as functions F^, F and Fc acting on non-amputated pomeron 
functions according to ( 1741) . To simplify notations we introduce a vector 



h(n,r,) = ^- /^ + % . (91) 
rf (ri + r2)^ 



From the results obtained in ^) we find 



rrf(ri|r2,r3) = -^^Vi 



|e«('-2-r3)h(^2, ri)h(r3, n) - e^"("^""^+''^^h(r2, ri)h(r3, -n) 
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1 - e'''''')—h(r2, ri)5^(r3 - + —(1 - e'''''')5\r2 - ri)5^(r3 - 

-V-'e'^''Vl6'{r2 - n6'{r, - n))}, (92) 



r(rih,r3) = ^Vt 



|e^«(r-2-r3)h(^2, ri)h(r3, n) - 2e"'^''-''+'-'Mr2: ri)h(r3, -n) 

Sime"""'^h(r2,ri)S\r3 - rs) - 3me'^''^(l - e'^'A^h{r2,ri)6\r3 - n - ra) 

-3e'^'^h{r2,n)h{r2 + rs,n) + 2me"'''[l - 2e'^'-^)^\i{r2,ri)5\r3 - n) 



+ — (1 + e''''^)6\r2 - n)8\r3 - n) - ^V-\'^'^Vl5\r2 - - n))}, (93) 

2asN^^ 
rc(ri|r2,r3) = -^2^Vi 

I _ ^i,^r,-r^+r,)^^^^^ n)\i{r3, -n) - ?,e'^''\v{r2, ri)h(r2 + n) 

-Ziire^^'^il - e"^^'^)4h(r2,ri)5(r3 - n - ra) + 2me''''^^\i{r2,n)5{r3 - rs) 

+2me^'"^2 4h(r-2,ri)5(r3 - n) - 167rVr'e^''"i5'(r2 - r^d'^ir^ - n) 

- ^V-^e-'-iV?52(r2 - ri52(r3 - r,)]. (94) 

Due to symmetry properties of the pomerons, these expressions have to be symmetrized 
in r2 and r^. Also all exponentials have to be substituted by their real and imaginary 
parts: 

e*'"'' cos kr, ie*'"' - sin kr, (95) 
which makes the cut vertexes real. 
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